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Abstract 

It is well known that high energy data alone do not discriminate between asymptotic In s and In 2 s 
behavior of pp and pp cross sections. By exploiting high quality low energy data, analyticity resolves this 
ambiguity in favor of cross sections that grow asymptotically as In 2 s. We here show that two methods 
for incorporating the low energy data into the high energy fits give numerically identical results and 
yield essentially identical tightly constrained values for the LHC cross section. The agreement can be 
understood as a new analyticity constraint derived as an extension of a Finite Energy Sum Rule. 



High precision low energy data represent a powerful constraint on the high energy behavior of hadronic 
cross sections via duality[Q[2|. The low energy data can be separated into two energy regimes, the resonance 
region and a region with energies in excess of a laboratory energy v$ where the resonances average into a 
featureless cross section in the sense of duality. These data represent powerful constraints on asymptotic fits 
to high energy data. Igi and IshidaT realized these constraints using a Finite Energy Sum Rule (FESR) 
which numerically averages the resonances, while Block and Halzenj2] simply required that the high en- 
ergy amplitudes fit both the experimental cross sections and their derivatives at the transition energy Uq. 
Both methods discriminate between a In s and In 2 s asymptotic behavior of the asymptotic cross section, 
conclusively favoring the latter. They appear to be more selective than conventional fitting techniques^. 

In this note we will show that the constraints of Block and Halzen[2] derive from analyticity^j, as does 
the FESR(2) of Igi and Ishida^Q. The purpose of this note is to show that they are in fact equivalent, as 
confirmed by fitting the two apparently very different methods to a common data set of pp and pp cross 
sections 0. 

Following Block and Cahn[7|, we describe the high energy data in terms of real analytic amplitudes 

/+ = 1 £{ A + PM*/*a) ~ W2] 2 + c^-V^ 1 -")/ 2 - »^/+(0) j , (1) 

for the crossing-even real analytic amplitude 8 , 9 and 

/_ = -^-D8 a - 1 e i7!{1 - a ^ 2 , (2) 
An 

for the crossing-odd real analytic amplitude. If a < 1, it parameterizes the Regge behavior of a crossing-odd 
amplitude which vanishes at high energies. A, a, (3, c, D, sq and \x are real constants. The variable s is 
the square of the center of mass system (cms) energy, v is the laboratory energy and p is the laboratory 
momentum. The real constant /+(0) is the subtraction constant required at v = 0. From the optical 
theorem we obtain the total cross section 



a ± =A + [3 



In s/ s - — 
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with p, the ratio of the real to the imaginary part of the forward scattering amplitude, given by 

p ± = — I f3n\ns/s - c cosinn^s^ 1 + —f + (0) ± Dsm(na/2)s a - 1 \ , (4) 
Otot L v ) 

where the upper(lower) sign refer to pp(pp) scattering. 

In the high energy limit, s — > 2mv where m is the proton mass, Eq. @ and Eq. Q can be written as 

o-of— ) = co + cilnf— ) +c 2 ln 2 (-) + fr, f—Y , (5) 
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where we have introduced the even high energy cross section cr [ — as weu as the relations A = c + — 
s Q = 2m 2 e~ Cl/(2c2 \ (3 = c 2 ,c = ^j-^j^Pv and D = f 7(J a /2) 5 - Tnis transformation linearizes Eq. © 
in the real coefficients Cq, ci, C2, /3p' and 6, convenient for a x 2 fit to the experimental total cross sections 
and p- values. Invoking Regge behavior, Igi and Ishida[Q fixed p = 0.5, which is the value we adopt in order 
to directly compare our method to their FESR(2) constraint. 

At the transition energy vq where we will match our high energy fits to the low energy data, we define 

&even(Vo) — 

= c + c 1 ln(Vm) + c 2 ln 2 (iVm) + /VK/m)' 1 - 1 , (8) 
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where a + (a ) are the total cross sections for pp(pp) scattering. Using the Block and Halzen valuepj] for 
fevenfVo), i-e., £r C vcn = 48.58 mb at vq = 7.59 GeV, we obtain the constraint 

co = crcvcn(^o) - ci ln(v /m) - c 2 lTL 2 (u /m) - /S-p^vo/m)^ 1 

= 48.58 - 2.091ci-4.371c2-0.3516/3p/. (9) 

In brief, we have used the pp and pp cross sections at the transition energy is — 7.59 GeV to anchor the 
asymptotic fit to the low energy data. The precise choice of v$ is not critical, as we will see further on. We 
will actually show further on that Eq. @ is a consequence of analyticity. 

To summarize, our strategy is to exploit the rich sample of low energy data just above the resonance 
region, but well below the energies where data are used in our high energy fit. At the transition energy vo, 
the experimental cross sections o~ pp (vo) and o~ PP (vo) are used to determine <7 eV en(^o) °f Eq. ©■ In turn, this 
constrains the asymptotic high energy fit so that it exactly matches the low energy data at the transition 
energy constraining the value of cq in Eq. @. Local fits are made to data in the vicinity of vq in order 
to evaluate the cross sections that are introduced in the above constraint equation, Eq. . We next impose 
the constraint Eq. || on a x 2 fit to Equations JfjJ) and (J7J). For safety, we start the data fitting at much 
higher energy, z/ m ; n = 18.72 GeV (^/Smjn = 6 GeV), well above vq. 

Given the previous analyses 012] we only consider an asymptotic In 2 s fit; the even amplitude parameter 
cq is constrained by Eq. (jSJ), i.e., by c±, c 2 and f3-p> and the experimental value of a even (po). We then perform 
a simultaneous fit to the experimental high energy values of a pp , a pp , p pp and p pp using six parameters: the 
even parameters c\, C2, (3-p> and /+(0) and the odd parameters S and a. Only the first 3 parameters are 
needed to describe the cross section. 

We now derive a new constraint from analyticity 4, , closely following Igi and Ishida'spP derivation of 
the fixed-energy sum rule 1 FESR(2) which they used to constrain their high energy fit. They wrote the 
imaginary part of their even high energy amplitude expressed in terms of their dimensionless coefficients 
Co, Ci, C*2, and B-p* as 

Im/+0) = ^- [Co + Cr In(>o/m) + C 2 ln 2 (v /m) + B v , (yo/m)"- 1 ] ■ (10) 
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After using the optical theorem, it yields the same cross section as oo of Eq. (JSJ) in the high energy limit 
(p — > v), if we make the substitutions Co — » ^-Co, c\ — ► t-Ci, C2 — > ^-C2, and [3-p> — > -^B-pi; note that 
the Igi and Ishida coefficients are dimensionless, whereas our coefficients Co, ci, C2, and ft-p* have dimensions 
of mb. We next introduce the true even amplitude F+(y) (which, of course, we do not know), along with 
the odd super-convergent difference amplitude 5 

!//+(!/) = v[F+(y) -U{v)\, (11) 

which, because of analyticity [S], satisfies the odd super-convergence relation 

vlmf + {v) dv = 0. (12) 

Note that vf+(v) satifies Eq. I|12fl . although neither the odd amplitude vF + (v) nor the odd amplitude vf+{v) 
alone necessarily satisfies it. Using Eq. (|12|) . we now write 



/>oo />oo 

/ v \vnF+{v) dv = / v Im/ + (v) dv. 
Jo Jo 



(13) 



Up until now, we have only used analyticity and the fact that the odd amplitude vf + {v) is super- 
convergent. However, super-convergence of vf + (v) also implies that f±(y), if it is a good representation of the 
high energy behavior, approaches sufficiently close to F + {v) at some transition energy i>o — taken as an energy 
somewhat above the resonance region — that the difference between these amplitudes is neglectable. Hence, 

1 We have changed their notation to conform with that used in the present paper, replacing F by /, the energy N by vq, and 
using capital letters for the real coefficients of their high energy parametrization , i.e., letting co — ► Co, ci — > G\, C2 — > C2 and 
/3-p/ — » B-pi . In what follows, m is the proton mass, p is the laboratory momentum and v is the laboratory energy. 
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the integrand of Eq. (|12[) is essentially zero for energies above Vq, allowing us to truncate the integration of 
Eq. 1|13[) at vq, thus obtaining the finite energy relation 



/ vlmF+(v)dv= / f + (v)dv = (14) 
Jo Jo 

for sufficiently large v . 

Applying the optical theorem to Eq. i|14|) . we have the relation 

I vlmF + (v) dv + — / vpa even (v)dv = / — - [Co + C\ \n{vjm) + C 2 ln 2 (^/m) 

JO 4?r Jm Jo m 

+B T ,(v /m)> J - 1 } dv. (15) 



Because vpa even (v) dv = p 2 cr even (p) dp, where po = V^o 2 — m 2 , this is the FESR(2) derived by Igi 
and Ishida which follows from analyticity, much as dispersion relations do. 

Again using the optical theorem, we now rewrite Eq. I|15(l in a more general form as 



/ vlmF + (v) dv + — / vpa even (v) dv = — / 
Jo 47r im 4tt J 



"a 
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v ao(v) dv. (16) 



Another constraint can be derived from the observation that above relation is satisfied for any vq in the 
energy region above the resonance region, where the cross section is smooth. In particular, Eq. (|16fl is valid 
at v + Avq, where Av is small, i.e, < Av <C v . Subtracting Eq. l(TH|l evaluated at v from Eq. ((TB|) 
evaluated at vq + Avq yields 

v +Av pv + Av 

pva even (v) dv = / v 2 a (v)dv, (17) 

VQ J V 

where cr even (v) is the value of the cross section at laboratory energy v and ctq(v) is the cross section at v 
obtained from the asymptotic amplitude. In the limit of Avq — > 0, Eq. (|17|) yields the analyticity constraint, 

CTevenM = — X (T (vq ) ~ Cr {v ) , (18) 
PO 

a relation good to ~ 0.04% forpo = 10 GeV, the value chosen by Igi and Ishida. Thus, we see that analyticity, 
which is the underlying fabric of the FESR(2) derived by Igi and Ishida 1 (Eq. (jlG^ . also requires that the 
high energy even cross section cto(^o) of Eq. (0 must match the low energy experimental cross section 
<Jeven{vo), if the parametrization of the high energy forward scattering amplitude is a good parametrization. 
Further, this result is independent of the value of the non-physical integral J Q vlvaF + (v) dv needed to 
evaluate the FESR(2) in Eq. (|16|l . even if it is very large. 

Igi and IshidaPP anchored their fit using Eq. H16[) . by numerically integrating p 2 a cvcn = p 2 \<Jp v -I- cr pp ]/2 
over the low energy resonance region below p$ — 10 GeV to obtain 



1 
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PO 

p 2 a even (p)dp = 3403 ±20 GeV. (19) 



Neglecting the error in Eq. H19|) and approximating the left-hand integral (the integral over the non- 
physical region) in Eq. (|15l) as 3.2 GeV, they obtained the constraint 

C = 8.87 - 2.04Ci - 4.26C 2 - 0.367B P /, (20) 

or, changing their dimcnsionlcss coefficients into our coefficients which have units of mb — by multiplying 
their coefficients by ^ — we rewrite the Igi and Ishida constraint as 

c = 49.28 - 2.04ci - 4.26c 2 - 0.367/V. (21) 

This is the constraint that will be used in an alternative fit to the high energy data. 

Before presenting our results, we comment on the "sieved" data that we will use for fitting[5]. It uses all of 
the data in the Particle Data Group 10 archive for pp and pp total cross sections and p-values with energies 
Js > 6 GeV. A robust In s fit was obtained which minimizes the Lorentzian squaredp], before imposing 
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the "Sieve" algorithm. The algorithm then proceeds iteratively to rid the data sample of the outliers, based 
on a maximum cut on the individual x 2 of the i th point, defined as A%|. Details are shown in ref. OEl- 
A value of x 2 /d.f.=5.657 was obtained for 209 degrees of freedom using the unscreened data^OI- This is to 
be compared to a value of x 2 /d.f.=1.095 for 184 degrees of freedom, when using a Axf max = 6 cut in the 
"Sieve" algorithmic . The "Sieve" algorithm eliminated 25 points with energies y/s > 6 GeV (5 a pp , 5 a pp , 
15 

Ppp): while changing the total renormalized x 2 from 1182.3 to 201.4. The 25 points that were screened out 
had a x 2 contribution of 980.9, an average value of 39.2. For a Gaussian distribution, about 3 points with 
Axf > 6 are expected, with a total x 2 contribution of slightly more than 18, not 980.9. This demonstrated 
the efficiency of the "Sieve" algorithm[B] in excluding outliers Ej. The same data set with Ax 2 max = 6 
and y/s > 6 GeV is used in the present analysis. 

Table ^ shows the results of a 6 parameter x 2 fit constrained by FESR(2) and, alternatively, by the 
analyticity constraint that matches <r cvcn at v§. The resulting x 2 have been renormalized [B] for the cut 
Ax 2 (x t ;a) = 6. Both fits are excellent, each with a renormalized x 2 per degree of freedom slightly less than 
1. 



Parameters 




a ~ In 2 s, A X 2 max = 6 






FESR2 Fit 


"Analyticity" Fit 


Even Amplitude 


c (mb) 




36.68 


36.95 


ci (mb) 




-1.293 ±0.151 


-1.350 ±0.152 


c 2 (mb) 




0.2751 ±0.0105 


0.2782 ±0.105 


/V (mb) 




37.10 


37.17 


M 




0.5 


0.5 


/(0) (mb GeV) 




-0.075 ±0.67 


-0.073 ±0.67 


Odd Amplitude 


S (mb) 




-24.67 ±0.97 


-24.42 ±0.96 


a 




0.451 ±0.0097 


0.453 ±0.0097 


X 2 ■ 




158.2 


157.4 


ft X Xmin 




180.3 


179.4 


degrees of freeedom (d.f) 




181 


181 


K x xLn/di. 


0.996 


0.992 



Table 1: The fitted results for a 6-parameter \ 2 fit with a ~ ln 2 (s) and the cut Ax? max = 6, for the FESR constraint 
c = 49.28-2.04ci -4.26c 2 -0.367,6V and the "analyticity" constraint c = 48.58-2.091ci -4.371c2-0.3516/3p/. The 
renormalized[S] Xmin/d-f-, taking into account the effects of the Ax? max cut, is given in the row labeled 1Z x Xmm/d-f- 
The errors in the fitted parameters have been multiplied by the appropriate r X 2 (see ref. |Sj). 

The pp and pp cross sections derived from the parameters of Table ^ are shown in Fig. j^i) as a function 
of the cms energy, y/s, for both methods. The pp (circles) and pp (squares) data shown are the "sieved" set. 
The short dashed and dot-dashed curves are the analyticity constraint fits to the pp and pp data, respectively. 
The solid curve and dotted curves are the some for the FESR fit. The difference between the two fits is 
negligible over the energy interval 4 < y/s < 20000 GeV; they agree to an accuracy of about 2 parts in 
1000. It should be emphasized that the FESR fit uses the experimental resonance data below y/s — 4 GeV 
for evaluating the constraint of Eq. (|21() , whereas the analyticity constraint fit uses the even cross section at 
y/s = 4 GeV for the evaluation of its constraint, Eq. (PJ, i.e., the alternative fits do not share any data. Both 
strongly support In 2 s fits that saturate the functional growth of the Froissart bound. 

In Fig. nj)) we show all of the pp and pp cross section data^U] m the cms energy interval 4 to 6 GeV, 
none of which was used in our high energy fit. Inspection of Fig. \T]p) reveals that we could have imposed the 
analyticity constraint anywhere from 4 GeV to 6 GeV without modifying the result. Thus, our conclusions 
do not depend on the choice of v$, the transition energy used in Eq. JHJl. 

Figure[21shows the fits for p pp and p pp as a function of the cms energy y/s; the "sieved" experimental data 
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are shown for y^s > 6 GeV. We conclude that the results are effectively the same for both fits and in good 
agreement with the experimental data. Accommodating p-values at lower energies allows one to constrain 
the cross section at higher energies by derivative dispersion relations, giving us additional confidence in our 
extrapolations. 

Summarizing, the FESR method and the new analyticity constraint introduced here yield fits to pp and pp 
cross sections and p-values that agree to 2 parts in f 000 over the large energy interval 4 GeV < y/s < 20000 
GeV. In particular, at the LHC energy of 14 TeV, the FESR fit predicts a pp = 107.2 ± f .4 GeV and 
p pp = 0.130 ± 0.002, whereas the analyticity fit predicts a pp = 107. 4± 1.5 GeV and p pp = 0.131 ± 0.002. We 
showed that this agreement was expected — it is numerical confirmation that analyticity, in its two guises, 
gives identical numerical results. Further, the fact that the renormalized \ 2 per degree of freedom in Tabled 
is excellent, giving a high probability fit, means that the choice of our high energy even asymptotic amplitude 
of Eq. (jSJ satisfies the analyticity constraint. It did not have to — had we used a poor representation for the 
even asymptotic amplitude, forcing the fit to go through the even cross section data at y/s = 4 GeV would 
have resulted in a very high x 2 - This was demonstrated in references [2] and pQ, where an asymptotic Ins 
paramctrization was decisively rejected. 

The fit of Block and Halzen[3] which additionally constrains the cross section differences, as well as 
derivatives of the cross sections at y/s — 4 GeV, for both pp and pp, yields essentially the same cross 
section and p-value, but with smaller errors. Clearly, from analyticity considerations, this technique is 
equivalent to evaluating additional FESRs, but is much more tractable numerically. This new tool yields 
both robust and precise values for the total cross section at the LHC energy of 14 TeV, as well as at cosmic 
ray energies, allowing us to make the prediction that at the LHC[2], p pp (14 TeV) = 0.132 ± 0.001 and 
a pp (U TeV) = 107.3 ± 1.2 mb. 
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Figure 1: a) The fitted total cross sections a vv and op v in mb, vs. y/s, in GeV, using the single constraint of Equations for 
the analyticity fit and J2D for the FESR fit of Tabled The circles are the sieved data for pp scattering and the squares are the 
sieved data for pp scattering for y/s > 6 GeV. The short dashed curve and dot-dashed curves are the analyticity fits — the even 
cross section at 4 GeV was fixed — to the pp and pp data, respectively. The solid curve and dotted curves are the FESR fits to 
the pp and pp data, respectively. It should be pointed out that the FESR and analyticity curves are essentially indistinguishable 
numerically for energies between 4 and 20000 GeV. b) An expanded energy scale that additionally shows the cross section 
data that exist HOl between 4 GeV, where <r cvcn was fixed, and 6 GeV, the beginning of the fitted data. It should be emphasized 
that none of the data between 4 and 6 GeV were used in the fits. We note that that the fits go through all of the unused points, 
with the exception of the pp point at 4.2 GeV, which would have been excluded by the "Sieve" algorithm 6 because of its large 
Ax?, had it been used. 
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Figure 2: The fitted p- values, p pp and p pp , vs. *Js, in GeV, using the single constraint of Equations J^J for the analyticity fit 
and 1211 for the FESR fit of Table HI The circles are the sieved data for pp scattering and the squares are the sieved data for 
pp scattering for y/s > 6 GeV. The short dashed curve and dot-dashed curves are the analyticity fits — the even cross section at 
4 GeV was fixed — to the pp and pp data, respectively. The solid curve and dotted curves are the FESR fits to the pp and pp 
data, respectively. It should be pointed out that the FESR and analyticity curves are essentially indistinguishable numerically 
for energies between 4 and 20000 GeV. 
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